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Flexoelectric deformation of a homogeneously polarized ball 
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The problem of finding the elastic strains arising due to the converse flexoelectric effect is solved 
for a case of uniformly polarized ball of isotropic dielectric material. The strains occurs due to the 
fact that, as shown in a previous paper, in presence of flexoelectricity elastic boundary condition 
has non-classical form. 

PACS numbers: 77.22.-d, 77.65.-j, 77.90. +k 



I. INTRODUCTION 



Direct flexoelectric eff'ect is the linear response of the electric polarization to a gradient of elastic strain. Although 
this effect has been known for a long time (see^ and references therein) , some of its features have been understood 
only recently. In particular, according to general principles of thermodynamics converse flexoelectric effect should 
also be present. At the same time, there is a certain asymmetry in the bulk constitutive electromechanical equations: 
while a linear strain obviously contributes to the homogeneous part of the polarization, homogeneous polarization 
ZAi ' does not appear in the equations that define the elastic stresses. Based on this, some authors have concluded that 
the flexoelectric effect, at least for a special arrangement, has no reversibility (see references in^). 

The above paradox is resolved as follows^i^. The equations of elastic equilibrium are differential, so they should be 
appended by boundary conditions. In the presence of flexoelectric effect elastic boundary conditions have non-classical 
. i form, and they include not only the gradient of polarization, but also the polarization itself^. This fact leads to the 
jrt ' homogeneous polarization distorting the body of finite size. 

In the practical sense the importance of flexoelectricity is that it, in contrast to piezoelectricity, is symmetrically 
allowed in centro-symmetric materials and therefore broadens the choice of materials which can be used for electrome- 
'^ I chanical devices. Besides, even though fiexoelectricity is a weak effect at the macroscale, at a nanoscale it becomes 
^ ' much stronger: reduced dimensions imply large gradients. This is why flexoelectricity has recently become a subject 
O ' of intensive studies. Particularly flexoelectricity in nanowires and nanopills was studied theoretically in^ but using 
classical elastic boundary conditions. 

Elastic boundary conditions obtained in^, were used in^i^, but only for the flat part of the body surface. At the 
same time, it was shown that the curvature of the surface plays an important role in the boundary conditions"^. It 
is interesting to solve the problem of flexoelectric deformation given a body surface curvature, which is the subject 
of this paper. Here we restrict ourselves to the case of the spherical shape of the body of an isotropic dielectric. 
Although this is a very special case, to solve the problem exactly for the more general case is difficult. The exact 
solution, except that it is interesting in itself, is useful as a test example in development of approximate methods of 
solution for more general conditions. 

Below, we also obtain differential equations and boundary conditions in a manner different from^, using curvilinear 
X^ , coordinates. It is easier than converting the previously obtained boundary conditions to these coordinates that appear 
most suitable for solving such problems. 



K^ ■ II. EQUATIONS OF EQUILIBRIUM AND BOUNDARY CONDITIONS IN CURVILINEAR 

Ut '. COORDINATES 

.9. 

Here we need a mathematical tool which is almost identical to the one used in general theory of relativity. The only 
difference is in the fact that, since the space is flat, the Riemann tensor is equal to zero, and the covariant derivatives 
commute. This mathematical tool is described in many standard textbooks, here we make only a few brief comments 
to introduce symbols and definitions. 

Curvilinear tensor indices are denoted by Greek letters. It will be convenient to assume that there is also a 
Cartesian coordinate system. Cartesian indices are denoted by Latin letters. For curvilinear indices is necessary 
to distinguish between top and bottom, for the Cartesian indices there is no need for that as they will always be 
bottom. For repeated indices, unless otherwise noted, the Einstein summation convention is adopted. The frame 
vectors are defined as Cia = dri/dx" = ri^a , where we use the notation of partial derivative by the index separated 
by a comma. Here r^ are the Cartesian coordinates of the point and x" are the curvilinear coordinates. Curvilinear 
tensor components are defined as Aa.../3 = Ai,,,jeia ■ ■ ■ Sjp . Raising of the Greek indices is performed in a standard 



way using the tensor g"" , which is reverse to the metric tensor g^p — eia&ifi ■ Lowering of the Greek indices is 
performed by gap ■ By g we denote the determinant of g^p . Invariant volume element is given by dV = y/gcfix , 
where d^x = dx^dx^dx^ . 

When equations in Cartesian components are transformed to equations in curvilinear components, the usual partial 
derivatives are replaced by covariant derivatives, which are denoted by lower Greek indices separated by semicolons. 
Covariant derivatives are defined as 

^ai...a„;/3 — ^ai...Q„,/3 ^7...a„-L ^iP ■ " ■ ^Qi...7-'- a^p j \^) 

where 

rI/3 = e]eia,p = -g'^^gscp + gsp^a - g^p.s) (2) 

are Christoffel symbols. 

Now we turn to the derivation of differential equations of equilibrium and boundary conditions. The free energy is 
written in the form 

F = [{Tel +Tp + Tficc)V9d^x , (3) 



where 



Tel = Iv'^^^'^'^Ua-p-eU^-S-X + \c''^^'ua-pU.,,S , (4) 

Tp = ]^a^^PaPp + \h^^'''Pc.;pPr..s " E^P^ , (5) 

^flx = \r^^\Po.-pUr,S - PaU^-5;p) . (6) 

The need for the higher terms in Tei which are proportional to Ua-.p-.e'i^-y-s-.i^ was shown ii*^. 

By varying the components of the polarization Pa with vanishing boundary terms, we obtain differential equations 
of polarization equilibrium: 

a^0p^ _ b'-MSp^,^,^ _ ^- „ r^^'u^;5-P = . (7) 

Similarly, varying the elastic displacements, we obtain differential equations of elastic equilibrium: 

'^"^;^=0, (8) 

where 



^Q/3 _ ^ap^S 



Uj-S 



■ r'"^Pr,5 - V'^^'^'^u^.s-ce ■ (9) 



To derive the boundary conditions for these differential equations when the surface is free, we require the vanishing 
of additional boundary terms arising from the variation. Boundary conditions for the equations of polarization 
equilibrium obtained immediately: 

(b''^'''Pr,snp + ^r^'''u^,snp^ =0, (10) 

where np is unit vector normal to the surface. 

As for the boundary conditions for the equations of elastic equilibrium, an integral equation quite similar to that 
derived in^ appear, and one needs to convert a surface integral. In^ such a conversion was done using additional 
surface coordinate system. In this case, when the curvilinear coordinates have been introduced from the beginning, 
any additional surface coordinate system is not necessary. It is sufficient to require that the equation of the body 
surface has the form: x^ = const . It is essential that the solution of specific problems is convenient to do in the 
coordinate systems of just this type. 



In this coordinate system, integration by parts is performed directly, and we obtain the elastic boundary conditions 
in the form: 

e"^'^n/jn^ = . (11) 

(^"7 _ 9-/37^^ + e-^^r]^ + (e-3^ V5) ,3 9-'^') % = , (12) 

where 

0a/37 ^ t;«'3e57C„^._5.^ _ If-f^t^Pg . (13) 

Naturally, these equations are only valid at the surface, we do not show it in the formulas for simplicity. 

The boundary condition ([TT|) is quite similar to that obtained in^. Otherwise the condition (|12p has a completely 
different look. Nevertheless, using PT|) and performing a series of transformations, (IT^ can be reduced to the form 
similar to the second elastic boundary condition in^: 

a^'^n^ - e°''^\j^n^ + e°'^''.sn^nfsn^ - e°'^'^-/fs^ = . (14) 

Here the tensor jp^ is expressed in terms of Christoffel symbols: 

7/3^ = r^^n'^n^ne + Tpgn^n.n-f - T^^n^ . (15) 

We emphasize that this relation is valid only in the coordinate system in which the equation of the surface is x^ — const . 
The fact that this equality cannot be valid in an arbitrary coordinate system is clear from the fact that the Christoffel 
symbols are not tensor components. Also it is possible to show that the equation (jlSp is consistent with the equation 
for 7q^ derived in^. Moreover, one can show that the tensor ja/s is symmetric, this was not obvious from the equations 
written in^. 

In general, equation ([7]) and (|S]) should be solved in the region occupied by the body with the boundary conditions 
([TOl) . (fTTj) and (fT4|) . In order to determine the electric field £'„ , the Poisson equation should be solved jointly, within 
the body as well as outside. Solutions of the Poisson equation inside and outside the body should be sewn together 
with the standard electrostatic boundary conditions. Electrostatic boundary conditions and Poisson equation in a 
curvilinear coordinate system are obvious, so we do not write them down. 

III. HOMOGENEOUSLY POLARIZED BALL OF ISOTROPIC DIELECTRIC 

Now we will apply what was said in the previous section to the case of an isotropic dielectric ball. We use the 
conventional spherical coordinate system: x = r sin 9 cos "tp , y — r sin sin ip , z = r cos . Here x , y and z are 
Cartesian coordinates, ip = x^ , = x'^ and r — x^ - curvilinear ones. Equation of the surface has the form r = R, so 
that these curvilinear coordinates belong to the special class described in the previous section. Instead of curvilinear 
indices 1 , 2 and 3 , we will use the indices ^p , and r , respectively. So the letters r , ip and are excluded from 
notation of "running" indexes. 

In Cartesian components the second rank material tensor of an isotropic medium reduces to a scalar: a,j — adij , 
where Sij is the Kronecker delta. Material tensors of the fourth rank are determined by two independent constants. 
So the elastic tensor can be written as Cijki — Ci25ij5^.i + CaidikSji + SuSjk) ^ the tensors 6^^; and fijki have a 
similar form. Sixth rank tensor of high elastic modulus of an isotropic medium is also determined by two independent 
constants: 

Vijlknm = Vi{5ij5lrn5nk + ^ij^nl^mk + SijSnmSlk + ^ikSjn^lm + SuSjnSkm + ^i7nSjnSkl + 

SinSjk^lm + SinSjlSkm + Sin5jm^kl)+ (16) 

V2{Sik^jlSnm + ^ikSjmSnl + SuSjkSnm + ^il^jm^nk + ^im^jk^nl + SimSjlSnk) ■ 

To write these tensors in curvilinear components one need only to replace the Kronecker delta Sij by the metric tensor 
with the corresponding location of the indices. 

Using the general formulas for calculating the components of the metric tensor, its determinant and the Christoffel 
symbols, and also using the material tensors of an isotropic medium, one can obtain partial differential equations of 
the problem by absolutely straightforward but cumbersome calculations. We will not write them down here because 
of their bulkiness and the fact that the calculations are quite standard. Instead, we will move to a discussion of 
variables separation. 



Generally the variable separation in such problems is performed by means of expansion of vectors in series in 
spherical vectors and expansion of scalars in series in scalar spherical harmonics. Note that, because of our use of 
covariant formalism here, the components of spherical vectors do not coincide with the ones usually used^. But the 
needed components can be easily derived from the conventional components. 

In our case the expansion in spherical functions is essentially simplified: a series in spherical harmonics contains 
only the terms with I = I , m = . Indeed, if far away from the ball the electric field is directed along the z axis and 
is homogeneous, then the electrostatic potential has only such a term. But then, since the variables are separated, 
the other variables also have only such terms. As a result, using the fact that lioCV'i ^) ~ cos 9 , it turns out that the 
angular dependence of the electrical potential is reduced to cos 6* ; the r-components have the same dependence, the 
^-components are proportional to sin0 and there are no ^/'-components. Thus, the problem is reduced to ordinary 
(radial) differential equations, and it is only necessary to provide a specified angular dependence explicitly. 

The radial system of equations and the boundary conditions for the complete radial problem (the ball in a uniform 
external electric field) is too cumbersome to write it here. Therefore, we restrict ourselves to a simplified version when 
the polarization is fixed and homogeneous. Thus only the elastic equations are solved. Since flexoelectric moduli of 
real materials are very small, this is a good approximation. 

In a simplified case, the system of radial differential equations has the form: 

vsi^fr + ^ef[" - 8ef[' + wi + ^er^ - i&h)+ 
v^ier + 4evr - e^vr + 12/1 - ^eiT + ^ef2 + 4^/2 - 12./2) = (17) 

C44e'(2eVr + ^U'l - 6/1 - 2U^ + &h)+ ci^ei^fi + 2^1 - 2/1 - 2^/^ + 2/2) , 



vsi^fr + ^en' - ^ef2 + 8/2 + ^eri - s/o - ^^4(2^/2' - 4/2 - efc - ^en + 2^/1 + 4/1 

C44?'(eV2' + 2^2 - 4/2 + U'l + 4/1) - Ci2e'(2/2 - U'l - 2/l) ■ 



(18) 



Here v^ = {vi + 2v2)R~^ , ^4 = (8wi -I- Av2)R^^ , prime denotes the derivative with respect to the dimensionless radial 
coordinate ^ = r/R, the components of the elastic displacement are expressed in terms of /i(^) and /2(0 as follows: 
Ur — /i cos 6 , ue — — -Rf/2 sin0 , m^ = . In the second equation we introduced a factor r = R^ so that /2(C) is an 
analytic function at zero for the final physical displacement in the center of the ball, as the g^^ = r~^ . The minus 
sign is added for convenience, so that /i = /2 = const means a homogeneous vector field directed along the z-axis. 
There are four boundary conditions for the radial equations: 



2^1 (/s" + I2 - W2 + 34/2 + 8/f + Uf[ - 34/1) + 2t;2(2/r + 2/^' - 20/^ + 86/2 
4/r + 16/( - 36/i) - 2c44i?'(/^ - /2 + h) = fi2R^P , 



(19) 



i?2[ci2(2/i - 2/2) + (ci2 + 2c44)/{] - vii9f[" + 18.f[' - 46/{ + 36/i - 14/^' + 32/^ - 36/2)- ,„„. 

2^2 (3/{" + 6/(' - 22/{ + 28/i - 2/^' + 16/^ - 28/2) = i?Vi2^ , ^ ' 

(18«i + 12«2)/{' -I- 12«i(3/{ - 4/1 - 2/^ + 4/2) = (/12 + 2/44)i?'P, (21) 

2vi{f!; + 2/^ - 6/2 + 2f[ + 6/1) + 4«2(/^' - 2/^ + 2/2 + 2f[ - 2/i) = U^R^P . (22) 

Note that the boundary conditions are non-uniform: the right hand side is proportional to the uniform polarization, 
which in this case is directed along the axis of z. The inhomogeneous boundary conditions leads to a deformation of 
the body in this case. 

IV. SOLUTION OF RADIAL EQUATIONS 

In general, for a system of two fourth-order equations (J17p - (J18p one needs eight boundary conditions, but here 
we have only four boundary conditions (J19p - (j22p . However, as usually in the case of such problems, the missing 
boundary conditions are replaced by the conditions that a solution should be analytic at zero. Therefore, the solution 
should be expressed in the form of series in non-negative powers of ^. However, it is more convenient to find the 
complete basis set of solutions, and then to present a general solution as a linear combination of such basis functions. 
If the number of basis functions is equal to the number of boundary conditions, then the coefficients in the linear 
combination are easily found by solving a system of linear algebraic equations. 



In accordance with the above, we express a solution in the form: 



/»(e) = E^fe^fe»(^)' 



(23) 



fc=i 



Bk^iO 



E 



(^kin<, 



(24) 



The constant coefficients akin are found from the condition that Bki {(,) obeys the system of differential equations p7)) 
- ([T5| . This condition leads to the following linear algebraic equations for the coefficients (index k which is numbering 
the solutions is omitted): 



[V3{n'^ - 2n^ - 9n^ + lOn + 16) + v4n'^ - 2n^ - In^ + 8n + 12)]ai„ + {vii^r? - 8n - 16) - v^{2n^' 
%r? + 2n + 12)]a2„ = [044(2^^ — 6n — 2) + Ci2(n^ — 3n)]ai„_2 — [c44(2n — 10) + Ci2(2n — 6)]a2„- 



(25) 



U3(4n^ - 4n - 8) + Vi{n^ + ri^ ~ An - 


-4)]ai„ + [v3(n4-2n3 


[C44(n + 2) + Ci2n]ai„_2 + [c44(»^^ - 


- 3n- 2) - ci22]a2„-2 • 



5n^ + 6r7, + 8) - Vi{2n'^ - 2n - 4)]a2, 



(26) 



From ([25|) - (|26|) one can see that the system of linear algebraic equations, even though it is infinite, has the 
characteristic structure: the coefficients with greater n are expressed in terms of the coefficients with smaller n. 
Therefore one can find ai„ for any finite n , one needs only to analyze the case of small n, and define some of the 
ain for these small n . The analysis shows that if one omits the physically meaningless solution corresponding to 
the translation of the ball as a whole along the axis z, then there are four independent solutions of (l^5|) - (I25|) i.e. 
four basis functions Bki{0 ; just the same number as the number of available boundary conditions. Note that one 
of the basis functions (say Bii{S,) ) can be chosen as a finite series that contains only the second powers of ^ . Such 
a function obeys not only the considered system of radial equations, but also the system of radial equations of the 
classical theory of elasticity. 




FIG. 1: Result of the calculation for i? = 1 • 10"* , P = 1 , C44 = 1.1 ■ 10^^ , C12 = 3.4 • 10^^ , /44 = /12 = 1 • 10"^ . Solid line is 
for t)i = 2 • 10"^ , -!;2 = 1 • 10"^ , dashed line is for ui = 2 • 10"^ , 112 = 1 ■ 10"^ . 



Without discussion of the following, rather conventional calculations, we give the numerical results (see the figure). 
Note that the graphs shows that, except for a thin layer near ^ = I , the curves are parabolic, the more so the 



less vi and t>2 are. In the bulk of the ball the only function Bii{(,) remains with a good approximation. Thus for 
a small Vi except for a thin surface layer, solution approximately obeys the equations of classical elasticity theory. 
This observation is very worth to note because it tells how one can built an approximate method for solving such 
problems, which uses the classical theory of elasticity in the bulk of the body. Certainly, in the thin layer near the 
surface anyway one should use equations of a theory of elasticity modified by flexoelectricity as described above. But 
a solution of this equations for a thin layer is simpler problem then the solution of them for whole body. There is no 
need for such approach for a ball, but it can be useful for more general problems. 

V. CONCLUSIONS 

In this paper we show how one can construct an exact, in the framework of the continuous medium theory, solution 
of the problem of flexoelectric deformation of homogeneously polarized ball of an isotropic dielectric. This specific 
example illustrates that the homogeneous polarization leads to flexoelectric deformation of the body, although the 
homogeneous part of the polarization does not enter in the differential equations of equilibrium. The reason for this 
deformation is related to non-classical form, in the presence of flexoelectricity, of the elastic boundary conditions, 
which was previously shown in^ in a general form. 

Even for such a simple geometry the solution is cumbersome and requires an introduction of non-standard special 
functions Bki{C) ■ However, the exact solution is useful, at least as a test example in the development of approximate 
methods and as a basis for the construction of a perturbation theory. 
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